Low-frequency Magneto-optical Spectra of Bilayer 

Bernal Graphene 



Y. H. Ho 1 ' 2 , Y. H. Chili 1 ** , D. H. Lin 2 , C. P. Chang 3 , and M. F. Lin 1 ^ 

department of Physics, National Cheng Kung University, Tainan, Taiwan 701 
2 Department of Physics, National Sun Yat-Sen University, Kaohsiung, Taiwan 804 
3 Ccnter for General Education, Tainan University of Technology, 710 Tainan, Taiwan 

August 6, 2009 

Abstract 

The low-frequency magneto-optical absorption spectra of bilayer Bernal graphene 
are studied within the tight-binding model and gradient approximation. The inter- 
layer interactions strongly affect the electronic properties of the Landau levels (LL's), 
and thus enrich the optical absorption spectra. According to the characteristics of 
the wave functions, the low-energy LL's can be divided into two groups. This divi- 
sion results in four kinds of optical absorption peaks with complex optical selection 
rules. Observing the experimental convergent absorption frequencies close to zero 
magnetic field might be useful and reliable in determining the values of several hop- 
ping integrals. The dependence of the optical absorption spectra on the field strength 
is investigated in detail, and the results differ considerably from those of monolayer 
graphene. 

PACS numbers: 78.20.Ls, 78.67.Pt, 73.20.At, 71.15.Dx 
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Since few-layer graphenes were produced by the mechanical friction 1 ' 2 and thermal 
decomposition, 3 ' 4 they have attracted considerable experimental and theoretical studies. 
Monolayer graphene (MG) and bilayer Bernal graphene (BBG; AB-stacked bilayer graphene), 
two of the simplest few-layer graphenes, display very different electronic features. MG 
exhibits linear energy dispersions near the Fermi level Ep = owing to the hexagonal 
structure. However, BBG possesses two pairs of parabolic bands 5 ' 6 owing to the interlayer 
interactions and stacking configurations. 7-9 These unusual electronic structures of few- 
layer graphenes have been experimentally verified through angular-resolved photoemission 
spectroscopy. 10 ' 11 

In the presence of a uniform perpendicular magnetic field B = B z, the magne- 
toelectronic properties of few-layer graphenes have been obtained by measurements of 
the quantum Hall conductivity, 2 ' 12 ' 13 scanning tunneling spectroscopy, 14 ' 15 and cyclotron 
resonance. 16,17 The theoretical results show that the linear bands of MG become the Lan- 
dau levels (LL's) characterized by E% v oc y/n c ' v B , 18,19 where c (v) is the index of the 
conduction (valence) LL's. E% v and n c,v are the energy and the integer quantum number 
of the nth energy band, respectively. The optical absorption peaks possess the selection 
rule | An| = \n c — n v \ = 1 and their frequency is proportional to \^Bo, 20 which has been 
confirmed by magnetotransmission measurements. 16 ' 21 ' 22 For BBG, the effective-mass ap- 
proximation, which neglects some interlayer interactions, predicts that the magnetic field 
changes the parabolic bands into LL's following E% v oc y / n c > v (n c ' v — 1)5 . 20 ' 23 ' 24 The opti- 
cal absorption frequencies obey a linear function of Bq in weak fields, and depend linearly on 



V^Bo as the corresponding energy leaves the parabolic band region. 20 The optical selection 

rule related to low-lying LL's is |An| = l, 23 the same rule that is applied to MG. How- 
ever, the calculations of Peierls tight-binding model including more important interlayer 
interactions indicate that the low-energy LL's should exhibit more complex features, 24 i.e., 
there exists an energy gap, asymmetry of LL's, two groups of LL's ( ls *LL's and 2nd LL's), 
and composite relations between E% v and B . The varied electronic features are reflected 
in the magneto-optical properties. 

In this work, the influence of the important interlayer interactions on the magneto- 
optical absorption spectra of BBG is investigated in detail by gradient approximation 25 ' 26 
and Peierls tight-binding model (see the geometry and detailed electronic properties in Ref. 
24). The band-like matrix 24 ' 27-31 is designed to solve a very huge Hamiltonian matrix. In 
this way, the wave functions 24 are obtained and can be used to divide the LL's of BBG into 
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two groups ls *LL's and 2nd LL's. For the magneto-optical excitations at a small magnetic 
field, the matrix will be larger and more complicated in the calculations, which makes 
the computation time very long. By utilizing the band-like matrix and obtaining the 
localized feature of wave functions, we can efficiently reduce the calculation time and thus 
evaluate the magneto-optical absorption spectra. The absorption peaks exhibit complex 
selection rules and composite field-dependent absorption frequencies. What deserves to be 
mentioned is that our study is more comprehensive than other studies in optical properties. 
The previous experimental and theoretical works mainly focused on the optical properties 
owing to ls *LL's. However, in our work, all of the possible optical absorption peaks resulting 
from the two groups are discussed in detail, and could be clearly distinguished from each 
other. Furthermore, a comparison of magneto-optical absorption spectra of MG and BBG 
is also made. 

The dispersionless LL's of BBG at B = 40 T are shown in Fig. 1(a). The interlayer 
interactions strongly influence the main features of LL's. The conduction and valence bands 
are asymmetric. There is an energy gap between the lowest unoccupied LL and the highest 
occupied LL. Based on the characteristics of wave functions, LL's can be classified into 
two groups, lst LL's (black curves) and 2rid LL's (red curves). ls *LL's appears at \E C ' V \ > 0, 
and 2nd LL's begins at E c > 7l + 7e (0.34 eV) and E v < - 7l + 7e (-0.39 eV). 7i 's 6 - 24 - 32 are 
the most important hopping integrals. The energies of these LL's exhibit the composite 



relationship of E$ v oc B Q and E$ v oc VBq 24 



The characteristics of wave functions deserve a closer investigation because they strongly 
dominate the optical properties. Each LL is composed of fourfold degenerate states with 
similar characteristics. At k x = k y = 0, the wave functions mainly consist of four localized 
subenvelope functions centered at positions x\ = 1/6, £2 = 1/3, ^3 = 2/3, and £4 = 5/6. 
%j=i-4 = A ijN /2R B or = B ijN /2R B , where i is the index of the ith layer and N represents 
the iVth Ai or Bi atoms. R 2 ^' 27 defines the dimensionality of the Hamiltonian matrix. The 
feature of wave functions resulting from the atoms with an even number (Af^ and B^) is 
similar to that with an odd {A ^ and B^) number. Because of this, only odd atoms are 
discussed in following calculations, where e (o) is an even (odd) integer. For convenience, 
one of the fourfold degenerate states at Xi = 1/6 is selected to be examined, as shown in 
Figs. l(b)-l(e). Through appropriate fitting, the wave functions of the nith ls *LL (black 
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curves) and n 2 th 2nd LL (red curves) can be expressed as 



A c l o , A c 2o oc <p 2 (x) , B c l o oc tp 1 (x) , B c 2o oc y? (x) , n{ eii = for m = 1; (la) 
^i,o> ^2,0 « V?o (a:) , B v l o oc y? 2 (x) , B v 2o oc </?i (x) , n\' eff = 1 for m = 1; (lb) 
A Z, A Z oc ip ni -i (x) , B C {1 oc ip ni - 2 (x) , B^ oc ip ni (x) for m ( = nf' elf ) > 2. (lc) 



A%, A% oc y? (x) , BZ oc y? 2 (x) , B c 2 % oc ^ (x) , n c /' e// = for n 2 = 1; (2a) 
^1,0, ^2,0 « <^n 2 -i (x) , 5^ oc v?„ 2 -2 (x) , B 2 f oc ^„ 2 (x) for n 2 ( = n 2 ,,; ' e// + 1) > 2. (2b) 

ip n (x) is the product of the nth-order Hermite polynomial and Gaussian function, 20 where 
n is the number of zero points in ip n (x) and chosen to define the quantum number of wave 
functions. 24 n°{ v ' e ^ (n 2 v ' e ^) is the effective quantum number of the nith lst LL's (n 2 th 
2nd LL's) defined by the quantum number of B 2 ^ ( A c {^ ). It should be noted that the 
ls *LL's and 2nd LL's could be distinguished by the different features of the wave functions 
in the mixed region (E c > 0.34 eV and E v < —0.39 eV), which is useful in analyzing the 
optical excitation channels. 

When BBG is subjected to an electromagnetic field at zero temperature, only excitations 
from the valence to the conduction bands occur. The optical selection rule Ak = due to 
the vertical transitions is mainly determined by the zero momentum of the photon. Based 
on the Fermi's golden rule, the optical absorption function is given by 
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where E is the unit vector of an electric polarization. The velocity matrix element M cv = 
(^ c (n^ 2 , k)|E • P/m e \ip v (n\ i2: k)) is calculated within the gradient approximation. M cv 
strongly depends on the characteristics of wave functions and would dominate the peak 
intensity. Due to the orthogonality of ip n (x), M cv has nonzero values except one condition 
when Af£ (Af^) and B\^ Q (B^) own the same <p n (x). Accordingly, the selection rules of 
optical absorption peaks can easily be obtained by the definition of wave functions in Eqs. 
(1) and (2). 

The low-frequency optical absorption spectrum of BBG at B = 40 T presents many 
prominent and some inconspicuous (indicated by arrows) peaks, as shown in Fig. 2(a). 



The prominent peaks with definite optical selection rules are discussed in the following 
paragraphs. The peaks can mainly be classified into four kinds of peaks, Wu's (black 
dots), u; 22 's (red dots), u^'s (green dots), and u^i's (blue dots), which originate in the 
excitations of ls *LI/s^ ls *LL c s, 2r ^LL , 's-> 2nd LL c s, ^LI/s^^LI/s, and 2nd LL^ ls *LL c s, 
respectively. The former (latter) two kinds of peaks result from the transitions between 
two LL's in the same (different) groups, as illustrated by the inset of Fig. 2(d). The 
four kinds of absorption peaks might have different optical selection rules because of the 
characteristics of wave functions in the two groups of LL's. 

The optical excitations of each prominent peak can clearly be identified, as shown in 
Fig. 2(a). For convenience, the excitations of wu's, u>22's, ^i 2 's, and u^i's are, respectively, 

1 X) £ J T £ £ T T D £ T T £ £ T T 1) £ T T £ £ T T 1 1) £ T T £ £ J T * 1 

represented as n 1 — > n{ , n 2 — > n 2 > , n{ — > n 2 , and n 2 — > n{ m the 
following. The selection rule of w^'s is denoted by Aria' (= nf, e ^ — n^ ,e ^). As to w u 's 
and u; 22 's, only the first peak (w^) of wu's originates in single channel, 1 — > 2. The other 
mth peak oj™ x (iv^i) consists of the pair oj™{ L and oj™{ H (c<%.' L and oo^ 11 ) corresponding to 
m m + 1 (m — 1 — > m) and m + 1 ^ m (m — > m — 1), respectively. The superscript L 
(if) indicates the lower (higher) energy of the pair. For example, lo\ x (^22) * s composed of 
o;^ 1, and cj^ r (u4 2 l and ^22) owing to 2 — > 3 and 3 — > 2 (0 — > 1 and 1 — > 0), respectively. 
The reason for this is that the energies of m — > m + 1 (m — 1 — > m) and m + 1 — > m 
(m — > m — 1) in w u 's (^ 22 's) are slightly different owing to the asymmetry of the conduction 
and valence LL's. In short, the optical selection rules of wu's and cj 2 2's can be represented 
by Ann = An 2 2 = ±1, which is the same for the LL's in MG. 

As for the transitions between two different groups, the mth peak of u^'s is formed 
with the pair, u^' 1 originating inm->m and uj^ H resulting from m + 1 — > m — 1. For 
instance, the first pair is indicated by lo\ 2 and uj\ 2 h in Fig. 2(a), which correspond to 1 — > 1 
and 2 — > 0, respectively. That is to say, cc^' L owns the optical selection rule An 12 = and 
ou 7 ^' 11 possesses Ani2 = —2. For u^i's, the excitations 0-^0 and — > 2 lead to the peaks 
UJ21 and of the first pair in Fig. 2(a), respectively. The other channels m — 1 — > m + 1 
and m — > m result in the pair c%[' L and lo™{ B \ respectively. In other words, the selection 
rules of U21S are An 2 i = and An 21 = 2. Obviously, the optical selection rules of u^s 
and U21S are different from those of w n 's and cj 2 2's, which is not clearly discussed in the 
previous theoretical works. 

The absorption rate significantly relies on the field strength and the excitation energy. 
In Figs. 2(a)-2(d), the peak intensity augments with increasing field strength, while the 
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opposite is true for the peak number. Furthermore, some peaks with commensurate ener- 
gies are overlapping because the peak spacing declines as the frequency increases, which 
varies the peak intensity. In contrast to BBG, MG exhibits absorption peaks with uniform 
intensity (Fig. 2(e)). The interlayer interactions are the main reasons for the difference 
between BBG and MG. The field-dependent absorption frequencies (a> a 's) of Wn's, w 22 's, 
cow's and a>2i's are shown in Fig. 3(a) by the black, red, green, and blue dots, respectively. 
cj a 's rise with increasing B and might follow the composite relationship of u a oc B and 
00 a oc v^o- Apparently, uj a oc y/B^ m MG (Fig. 3(b)) is very different from that in BBG. 
In addition, the convergent frequencies of Wn's, u; 22 's, u^'s and cj 21 's at B — > are ap- 
proximately 0, 0.73 eV (271), 0.34 eV (71 + 7 6 ), and 0.39 eV (71 — 7 6 ), respectively. This 
implies that the optical measurements can reasonably determine the values of 71 and 76 
through observing the convergent frequencies of absorption peaks. The predicted result is 
very useful and reliable for experimental studies. 

The tight-binding model is widely applied to tackle the physical problems of semicon- 
ductors and carbon-related systems. Through the band-like matrix and the localized feature 
of wave functions, the huge and complicated Hamiltonian matrix can be further simplified 
in the calculations of the magneto-optical excitations for small field cases. Thus the com- 
putation time is substantially reduced. Moreover, the numerical strategy can be utilized 
to study not only the electronic and optical properties, but also the Coulomb excitations 
and transport properties in other crystal materials. 

In summary, it can be said that the interlayer interactions strongly affect the magneto- 
optical properties of BBG in the presence of uniform perpendicular magnetic fields. Based 
on the characteristics of wave functions, the conduction and valence LL's are asymmet- 
ric and can be divided into ls *LL's and 2rid LL's leading to four kinds of optical absorption 
peaks Wn's, cj 22 's, u^'s, and 0^2 's. Their optical selection rules are An n = ±1, An 22 = ±1, 
An 12 = 0, —2, and An 2 i = 0, 2, respectively, which can again be explained by the char- 
acteristics of wave functions. The selection rules and field-dependent absorptions of BBG 
are more complicated than those of MG. The above-mentioned influences of interlayer in- 
teractions on the magneto-optical properties could be confirmed by magnetoabsorption 
spectroscopy measurements. 17,21,22 Furthermore, the convergent absorption frequencies at 
nearly zero field might be helpful and reliable in determining the interlayer atomic inter- 
actions 71 and 7 6 through the magneto-optical measurements. 
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Figure Captions 

FIG. 1. (a) Landau levels of bilayer Bernal graphene at B =A0 T. The wave functions 
of (b) A lj0 , (c) Bi , (d) A 2>0 , (e) B 2fi atoms with odd integer indices are shown. 

FIG. 2. The optical absorption spectra of bilayer Bernal graphene at (a) 40 T, (b) 30 
T, (c) 20 T, and (d) 10 T. The spectrum of monolayer graphene at 40 T is plotted in (e). 

FIG. 3. The field-dependent optical absorption frequencies of (a) bilayer Bernal graphene 
and (b) of monolayer graphene. 
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